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Abstract 

Alfven's  perturbation  method  for  treating  the  motion  of  a  point 
charge  in  a  magnetic  field  consists  of  three  elements:  l)  The  field  is 
regarded,  in  zero  order,  as  homogeneous.  In  the  next  order,  the  zero- 
order  motion  is  replaced  by  its  averaged  effect,  a  magnetic  dipole,  and 
the  interaction  of  this  'equivalent  dipole*  with  the  inhomogeneities  of 
the  field  constitvrtes  the  perturbation,  2)  The  moment  of  the  'equivalent 
dipole'  is  constant*  3)  If  the  field  is  axially  symmetric,  the  'equivalent 
dipole'  moves  on  the  surface  generated  by  rotating  a  line  of  force  about 
the  polar  axis. 

Alfven  develops  his  method  by  the  analysis  of  several  special  cases. 
In  this  report,  we  have  re-derived  his  perturbation  scheme  from  a  viewpoint 
that  emphasizes  general  dynamical  principles.   Thus:  1)  is  derived  directly 
from  Hamilton's  principle,  2)  is  proved  by  the  adiabatic  theorem,  after  not- 
ing that  the  moment  of  the  'equivalent  dipole'  is  an  action  vsiriable,  while 
3)  is  related  to  the  conservation  of  angular  momentum. 

Our  remarks  under  2)  and  3)  are  intended  as  an  extension  of  Alfven's 
discussion. 
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1,  Introduction 

The  basic  elementary  process  in  magnet o-hjrdrodynamics  is  the  motion 
of  a  point  charge  in  a  magnetic  field.  For  a  uniform  magnetic  field,  the 
solution  is  trivial;  for  a  more  general  magnetic  field,  the  solution  can  be 
quite  complicated.  An  intermediate  situation  has  been  considered  by  Alfven*-  ^ 
namely,  the  case  where  the  variation  of  the  field  may  be  regarded  as  a  small 
perturbation  of  a  uniform  field. 

In  Alfven's  method,  the  field  is  regarded,  in  zero  order,  as  homogen- 
eous. In  the  next  order,  the  zero-order  motion  is  replaced  by  its  average 
effect,  a  magnetic  dipole,  and  the  interaction  of  this  'equivalent  dipole' 
with  the  inhoraogeneities  of  the  field  constitutes  the  perturbation. 

Alfven  presents  his  method  by  analyzing  several  special  cases.  In 
this  report,  we  have  re-derived  his  perturbation  scheme  from  a  viewpoint  that 
emphasizes  general  dynamical  principles. 

Our  discussion  extends  Alfven's  treatment  in  two  respects:  First, 
the  circumstances  imder  which  the  moment  of  'equivalent  dipole'  is  constant 
are  summarized  concisely  by  noting  that  it  is  an  action  variable,  and  there- 
fore is  unaffected  by  any  changes  that  are  'adiabatically  slow'  compared  to 
the  zero-order  motion.  Second,  if  the  field  is  axially  symmetric,  it  is  shown 
that  conservation  of  angular  momentum  requires  that  the  'equivalent  dipole' 
must  move  on  the  surface  generated  by  rotating  a  line  of  a  force  about  the 
polar  axis. 

2.  The  zero-order  motion 

We  consider  a  point  particle  (charge  e,  mass  m)  in  a  magnetic  field,  H. 
The  position  vector  of  the  particle  is  the  sum  of  a  slowly  varying  part,  R,  and 
a  rapidly  varying  part,  r.   r  is  used  to  describe  the  zero-order  motion  and  R 
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to  describe  the  first-order  motion.  In  an  homogeneous  field,  R  is  constant. 

The  zero-order  motion  regards  the  magnetic  field  as  approximately 
homogeneous.  The  particle  then  makes  excursions  of  the  order  r  about  a  point 
R,  Thus,  the  approximation  of  homogeneity  is  justified  if 


(2.1)  L  •  'r  ^^W  «  S.(S.) 


If 


H(R)<>'-i- 
—     n 

r" 


(2»1)  requires  that 


i  »n 

r 

or,  since  n  ^  1  in  most  cases,  our  perturbation  theory  will  be  applicable  if 

(2.2)  I  »  1  . 

So  far,  we  have  held  R  constant  during  the  motion  and  have  considered 
the  variation  of  the  field  as  r  varies.  Actually  R  also  varies  if  the  field 
is  not  homogeneous.  If  the  consequent  magnetic  field  variation  is  to  be  negli- 
gible, it  is  sufficient  to  add  to  the  condition  (2.1): 

(2.3)  T  S.*^  L 

where  t'  is  the  period  associated  with  the  zero-order  motion.  (2.3)  implies 
that 

(2.U)  R  «  f 
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while  (2.2)  and  (2.3)  imply  that 

(2.5)  T2.  «  5.  . 

These  conditions  define  the  domain  of  applicability  of  the  perturbation 
method.  From  (2.U)  and  (2.5)  we  see  that  we  may  consider  R  ^   const,  during 
times  of  the  order  of  T'  • 

We  now  exhibit  the  consistency  of  these  restrictive  conditions  by  show- 
ing that  they  do  indeed  lead  to  the  type  of  zero-order  motion  we  wish  to  obtain. 

For  this  purpose,  we  start  with  the  Lagrangian 

(2.6)  L  -  ^  (r  +  1.)^  +  f  (r  +  R)  •  A(r  +  R)  . 

We  have  seen  that  R  ^   const,  during  times  of  the  order  of  'Y  •     Therefore,  if 
we  limit  ourselves  (in  Hamilton's  principle)  to  variations  of  the  trajectory 
that  are  different  from  zero  only  dviring  times  of  the  order  of  ^  ,  we  may 
replace  (2.6)  by 

(2.7)  L-'l  r^  +  f  £•  A(r  +R) 

where  only  r  is  varied,   and,   after  the  variation,  we  use   (2.2). 

Thus,   using  (2.7)  in  Hamilton's  principle  and  varying  r  only,  we  find: 


(2.6)  -f  -    -^    f   X  H(r  +  R) 

—         mc      —        —  — 


which,   according  to   (2.2),   is  equivalent  to 
(2.9)  r-    ^    tx-H(R) 

where 

n  (R)  -  -  -saiBi    . 

- —  ^—  mc 
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These  equations  are  solved  by: 

(2.1C)  £  =  X1(5.)>'£    . 

That  is,  the  particle  revolves  in  a  circle,  of  radius  r,  about  the  field  lines. 
The  angular  velocity  is  44 . 

Thus,  the  restrictive  conditions  we  have  postulated  do  lead  to  a  zero- 
order  motion  that  takes  place  in  a  homogeneous  field.  We  note,  from  (2.10), 
that  the  zero-order  motion  is  entirely  transverse  to  the  magnetic  field.  In 
zero  order,  then,  the  motion  along  the  magnetic  field  is  undetermined.  We  shall 
show  below  that  the  first-order  equations  of  motion  fix  the  motion  along  the 
magnetic  field  lines  completely. 

We  return  to  the  Lagrangian  (2.6).  This  time,  we  wish  to  vary  R.  Accord- 
ing to  (2.5)  only  variations  of  R  that  involve  many  periods,  '^  ,  of  the  zero- 
order  motion  will  be  noticeable.  That  is,  during  a  variation  of  R^  we  may  re- 
place quantities  involving  r  or  f  by  their  average  values.  We  therefore  pro- 
ceed to  the  calculation  of  the  required  averages. 

First,  we  note  that  the  term  n^  /2  in  the  Lagrangian  may  be  ignored, 
since  only  terms  involving  R  and  R  will  be  varied.  Next,  we  note  that  (2,1) 
implies  a  similar  restriction  on  the  vector  potential: 

r  .  7j^  A(R)  «  A(R)  . 

Finally,  since  the  zero-order  motion  given  by  (2,10)  is  such  that  the  particle 
executes  a  circle  about  the  field  lines,  we  have: 

r  "  0  »  ?  ♦ 

From  these  remarks,  the  effective  Lagrangian  for  variations  of  R  is 

(2,11)         L  --5  R^  *  f  S.  •  W  ♦  f  L  •  'r  A(R)  •  r 
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where  the  bar  indicates  an  average  over  the  zero-order  motion, 

3»  The  'equivalent  dipole' 

We  shall  now  show  that 


(3.1)  f  I  '  ^R  i^^^  •  £  ■  ti  •  11(2.) 

where 

E  -  27  ^=Q^(^) 

is  the  magnetic  moment  of  the   'eqtiivalent  dipole'. 

For  the  left  side  of   (3.1),  we  have,   using  (2.10), 

^TL^   \  A(R)  -  I     XLxTT  :   7p  A(R)    . 

We  therefore  have  to  compute  the  dyadic  r  r^  where  r  is  a  vector  in  the  plane 
perpendicular  to  the  field  H^  and  the  average  is  taken  over  a  circle  in  this 
plane.  We  shall  show  that 

(3.2)  n  '  T  h 

where  I.    is  the  unit  dyadic  in  the  plane  perpendicular  to  H. 

To  prove  (3.2),   we  note  that  r  r  is,   by  definition,   the  two-dimensional 
dyadic  r  r  averaged  over  all  directions  in  the  plane.     The  only  dyadic  that  is 
independent  of  direction  in  the  plane  is  I.j  therefore  r  r  is  proportional  to  I.. 
The  constant  of  proportionality  is  evaluated  by  requiring  the  spur  of  both  sides 
of   (3.2)  to  be  the  same.     This  proves   (3.2). 

On  the  other  hand,  if  a  cxirrent  I  threads  a  closed  contour  C,  the  mag- 
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2     I     £  ^  ^ 


In  our  problem. 


I-  ^ 
ct 


Hence,   iising   (2.10), 

r 

(3.3)  ii"^'^    r   Lxrdt 

e    .  "2 


^.r'^II(R)    . 


Assembling  these  results,  we  see  that  the  left  side  of  (3»l)  leads  to  the 
right  side. 

Returning  to  (2.11),  the  effective  Lagrangian  for  variations  of  R 
becomes: 

(3.U)  L  ^  I  r2  +  I  R  .  A(R)  +  M   o   H(R) 

-  ™  r2  +  £  R  .  A(R)  -  n(R)  H(R)   . 

In  the  s  econd  form,  where  the  scalars  |j.  and  H  appear,  we  have  used  the  fact 
that  ^  and  H  are  anti-parallel, 

U»  The  constancy  of  h.(R) 

According  to  (3.3),  the  magnetic  moment  of  the  'equivalent  dipole'  is 
a  function  of  R.  We  would  therefore  expect  p,(R)  to  give  a  contribution  to 
(3.U)  when  R  is  varied.  In  actuality,  however,  this  is  not  the  case:  p.  may  be 
treated  as  a  constant  in  (3.U)  if  the  variations  of  R  are  siifficiently  slow. 
We  shall  now  carry  out  the  reasoning  that  leads  to  this  conclusion. 
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We  first  consider  a  homogeneous  field  along  the  z-axLs,  ajid  we  intro- 
duce the  cylindrical  coordinates  (p,©, z).  The  z-coraponent  of  (canonical) 
angular  moraenttan,  P_,  is  constant.  That  is, 

(U.l)  Pq  -=  mp2©  +  i  Agp 


■  mp  6  +  s—  Hp  =  const. 
But,  from  (2,9), 


mc 


Therefore, 

(U.2)  p^  »  -  2|!.  JI  .  -  I^  ^,  -  const. 


The  action  variable  associated  with  the  coordinate  9  is 


2n 
>o 


(U.3)        Jq  -   Pe  '^^  "  2"  ^« 


=  -^2n,  . 

For  a  homogeneous  field,  therefore,  the  magnetic  moment  of  the  equivalent 
dipole  is  a  constant  of  the  motion  and  is  an  action  variable.  It  follows  from 
the  adiabatic  theorem  that  it  is  also  an  adiabatic  invariant,  that  is,  the 
magnetic  moment  is  unaffected  by  all  changes  that  are  'adiabatically  slow* 
compared  with  the  period  of  the  motion  corresponding  to  an  homogeneous  field. 

According  to  (2.5^),  we  have  limited  ourselves  to  motions  of  R  that  are 
negligibly  small  in  a  period  of  the  zero-order  motion.  Thus,  the  possible 
variations  of  R  in  (3.h)  will  satisfy  the  requirements  of  the  adiabatic  theoran. 
For  our  applications,  therefore,  \i  may  be  considered  a  constant  when  R  is  varied 
in  (3.U). 
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For  completeness,  we  note  that  ^  is  related  to  the  flux  encircled 
by  the  particle,  f,  by: 

Thus,  the  flux  encircled  by  the  particle  during  a  period  of  the  zero-order 
motion  is  also  an  adiabatic  invariant. 

As  a  final  remark  in  this  section,  we  note  that  VT  ,  the  mean  energy 
of  motion  transverse  to  the  field,  is  given  by: 

(U.6)  Wj^"-ti»H-^r  SL    "^r 

On  the  other  hand,  the  total  energy,  E,  is  given  by: 


That  is,  motion  along  the  lines  of  force  is  associated  with  R.  Although  this 
already  implies  that  the  'equivalent  dipole'  moves  so  as  always  to  remain  on 
the  same  line  of  force,  we  shall  prove  this  more  directly  in  the  next  section. 

5.  The  orbit 

We  introduce  a  system  of  spherical  coordinates,  (R,©,J^)  and  we  suppose 
the  vector  potential  to  have  only  the  component  A^  different  from  zero.  Since 
the  vector  potential  is  chosen  divergenceless,  it  is  independent  of  J^,  and  we 
are  considering  a  situation  where  the  field  is  unaffected  by  a  rotation  about 
the  polar  axis*  It  follows  that  Pj,  the  (canonical)  angular  momentum  about  the 
polar  axis,  is  a  constant  of  the  motion. 

Introducing  spherical  coordinates  into  (3.U),  we  find: 

(5.1)  P^  ■  II  «  m  R  sine  <R  sinO  i  *  ^  Aw>  -  const. 
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Since 

the  second  term  in  ($•!)  is  of  the  order 

22  R  -  nr(S.)  -  H-)  . 
mc         ^r     r 

The  first  terra  in  ($•!),  however,  is  of  the  order  of  R«  Therefore,  according 
to  (2.2)  and  (2.U)»the  second  term  is  very  much  larger  than  the  first  and 
(5»1)  becomes: 

($.2)  ^0  '^  I  V  ^  ^^"®  "  °°^^*» 

The  'equivalent  dipole'  therefore  moves  so  as  to  satisfy  (5«2)« 

To  see  ^riiat  this  means  geometrically,  we  observe  that  the  field  has  the 
components 


(5-3)         «<,'  -ffSS-A'V'^' 


"« 


0  . 


If  we  make  a  displacement  along  a  line  of  force,  the  components  5R  and  6©  of 
the  displacement  are  related  to  the  field  ty: 

5R     "r 
rS5    ^      * 

Substituting  from  (5»3)  in  this  relation,  we  findx 
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Thus,  (5»2)  holds  if  and  only  if  the  'equivalent  diple'  always  moves  on  the 
surface  generated  by  rotating  a  line  of  force  about  the  polar  axis. 
In  spherical  coordinates,  (3»U)  has  the  form: 

(5.5)  L  -  I  (ft2+  r2  ©2)  +  m  R  sine  0  k^  ^  *  jfe  A^>  '  ^    ' 

This  Lagrangian  corresponds  to  a  particle-charge  e,  mass  m,  magnetic  moment 
|j,  -  moving  in  the  magnetic  field  H. 

By  the  same  reasoning  that  led  from  (5«l)  to  (5.2),  (5»5)  may  be 
replaced  by 

(5.6)  L  »  I  (R^  +  R^  ©^)  +  I  A^  R  sin©  0  -  nH  . 
The  energy  integral  corresponding  to  (5»6)  is*. 


(5.7)  E-|  (r2  +r2  9^)+  nH  . 


The  actual  orbit  of  the  'equivalent  dipole'  may  be  obtained  as  follows: 
From  (5.6)  we  obtain  an  equation  of  motion  for  each  of  the  three  coordinates. 
The  equation  for  f6   integrates  at  once  to  give  (5.2).  As  already  discussed,  this 
prescribes  the  surface  on  which  the  'equivalent  dipole'  moves.  The  form  of  this 
surface,  as  given  by  (5«2),  may  be  written 

(5.8)  R  -  f(©). 

Putting  this  into  (5»7),  results  in  an  equation  for  ©.  Finally,  using  this  new 
equation  to  eliminate  6  and  §  from  the  equations  of  motion  for  the  R  and  ©  co- 
ordinates,  previously  derived  frcm  (5.6^  we  obtain  an  equation  for  0,     Integra- 
ting  these  equations  for  ©  and  0,   and  combining  with  (5.8),  specifies  the  orbit 
completely. 
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